from random import uniform

def cont_integral(a, b, N):
	summ = 0.0
	summ_n = 0.0
	for i in range(N):
		x = uniform(0, 1)**(1./(a+1.))
		area = x**(b-a)
		summ += area
	return summ/float(N)

a = 2.
b = 3.
N = 1000

print cont_integral(a, b, N)
